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Introduction 

Tumor  progression  depends  on  the  ability  to  stimulate  the  development  of 
new  blood  vessels.'  This  process  is  known  as  tumor  angiogenesis.  It  has  been 
hypothesized  that  the  characteristic  new  vessel  growth  associated  with  tumor 
angiogenesis  should  allow  differentiation  between  malignant  and  benign  breast 
lesions.  Recent  studies  point  to  vascular  endothelial  growth  factor  (VEGF)  as  a  major 
regulator  of  physiological  angiogenesis.  VEGF  is  overproduced  by  tumor  cells  and 
acts  predominantly  as  a  potent  angiogenic  factor  contributing  to  malignant 
progression.^  The  VEGF  receptors  are  specifically  expressed  in  the  cell  surface  of 
vascular  endothelial  cells.  Benign  and  malignant  tumors  are  characterized  by 
differential  cellular  expression  of  VEGF.  Monitoring  the  expression  of  VEGF 
receptors  may  therefore  be  a  reliable  marker  of  the  proliferative  status  of  tumors. 

Attempts  to  improve  tumor  detection  by  non-radioactive  magnetic  resonance 
imaging  (MRI)  techniques  have  led  several  investigators  to  develop  “immuno-specific” 
paramagnetic  contrast  agents.  We  propose  that  the  use  of  a  magnetic  label  with  high 
affinity  and  selectivity  for  VEGF  receptors  will  allow  VEGF-induced  microvascularity  to 
be  visualized  non-invasively  using  MRI. 

The  concept  of  this  research  is  to  develop  a  non-invasive  observation  of 
angiogenesis  through  the  design  and  validation  of  MRI  pulse  sequences  and  novel 
contrast  agents  and  to  use  this  agent  to  monitor  angiogenesis  in  vivo  over  time 
during  tumor  development. 
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The  specific  objectives  are: 

1 .  To  develop  novel,  highly  specific  iron  oxide  based,  magnetic  resonance  contrast 
agents  to  be  used  for  targeting  the  macrophage  population  associated  with  tumor 
cells. 

2.  To  optimize  MR  pulse  sequence  technology  for  the  relative  assessment  of  tissue 
iron  concentrations  based  on  the  quantitative  evaluation  of  T2  and  T2*  weighted 
images. 
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Body 

Characterization  of  Breast  Tumor  Animal  Model  (Phase  I) 

The  first  phase  of  this  research  involved  characterizing  the  histologic  and  MRI 
appearance  of  mammary  tumors  in  the  chosen  animal  model. 

Breast  tumors  were  induced  in  30-day-old  female  Sprague  Dawley  rats  by 
intraperitoneal  (i.p.)  injection  of  the  carcinogen  N-ethyl-N-nitrosurea  (ENU:  Isopack, 
Sigma  Chemicals,  St.  Louis,  MO).  ENU  is  known  to  be  a  potent  carcinogen  that  induces 
well  vascularised  mammary  tumors  in  rats.^  This  tumor  model  was  chosen  because  a 
range  of  benign  and  malignant  tumors  develop  depending  on  the  dose  of  ENU 
administered  and  because  the  tumor  morphology  and  angiogenic  behavior  closely  mimic 
breast  tumors  in  humans.  ENU  was  injected  at  two  different  doses:  40  mg/kg  (n=20)  and 
180  mg/kg  (n=20).  Animals  were  visually  checked  for  tumor  growth  ever  second  day. 
Once  tumors  were  identified  their  size  was  tracked  weekly  using  calipers.  Tumors  were 
imaged  once  they  reached  a  size  of  at  least  4  cm.  Animals  were  anesthetized  for  MRI 
with  ketamine/xylazine  (10/100  mg/kg,  i.p.).  High-resolution  in  vivo  MR  images  were 
obtained  of  all  mammary  tumors.  After  imaging  animals  were  sacrificed  and  the  tumors 
removed  for  histopathologic  analysis.  Photographs  of  the  position  of  the  tumors  in  the 
mammary  tissue  were  obtained  at  post  mortem  to  guide  the  correlating  of  MR  and 
histologic  images. 

Tumors  appeared  as  early  as  60  days  post  ENU  administration.  In  the  low  dose 
group  (40  mg/kg)  only  6/20  rats  developed  mammary  tumors  by  100  days  post  ENU 
administration.  All  6  of  these  tumors  were  classified  as  benign.  In  the  high  dose  group 
(180  mg/kg)  19/20  rats  developed  mammary  tumors;  12  animals  developed  multiple 
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mammary  tumors,  in  total  44  tumors  were  removed  and  examined.  Of  the  44  tumors 
from  this  high  dose  group,  35  were  classified  as  malignant  and  9  were  classified  as 
benign.  Based  on  the  results  of  this  experiment  further  studies  used  only  the  high  dose 
(180  mg/kg)  of  ENU  since  this  protocol  provided  a  useful  spectrum  of  tumor  types  for 
analysis. 

The  appearance  of  benign  mammary  tumors  in  the  MR  images  was  typically  of  a 
consistent  nature.  Benign  tumors  appeared  to  have  homogeneous  signal  intensity 
throughout.  Tumor  nodules  were  clearly  visualized.  The  MRI  appearance  of  malignant 
tumors  v/as  more  variable.  Malignant  tumors  often  appeared  to  have  a  localized  region 
of  very  low  signal  intensity  that  represented  a  necrotic  core.  Based  on  the  results  of  this 
analysis  the  tumor  type  and  regions  of  tumor  anatomy  could  be  recognized  in  vivo  by 
MRI. 


Pulse  Sequence  Development  (Phase  ID 

The  second  phase  of  this  project  involved  MRI  pulse  sequence  and  software 
analysis  development.  This  work  was  critical  to  the  concept  of  quantifying  iron  content 
in  vivo  using  noninvasive  MRI  techniques.  This  phase  of  research  was  done  in 
collaboration  with  Dr.  Brian  Rutt’s  group,  a  collaborator  in  the  Imaging  Research  Lab. 

3D  multi-echo  gradient  echo  and  spin  echo  (CPMG)  sequences  have  been 
developed  and  implemented  in  our  lab  (detailed  descriptions  of  the  pulse  sequence 
development  and  testing  can  be  found  in  the  appended  manuscript  (“Application  of  Static 
Dephasing  Regime  Theory  to  Superparamagnetic  Iron-Oxide  Loaded  Cells”).  These 
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experiments  have  resulted  in  the  submission  and  acceptance  of  a  manuscript  that  is 
currently  being  revised  for  publication  (see  Appendix). 

In  Vitro  Cell  Studies  (Phase  III) 

The  third  phase  of  this  work  involved  an  examination  of  iron  uptake  by 
phagocytic  cells.  In  invasive  breast  cancers,  the  neoplastic  cell  population  is  often 
outnumbered  by  stromal  cells,  such  as  the  tumor  associated  macrophage  (TAM),'^  a 
mononuclear  cell  of  the  inflammatory  system.  It  is  thought  that  monocytes  in  the 
peripheral  circulation  are  actively  recruited  to  the  tumor  site.  Once  recruited,  the 
monocytes  differentiate  to  become  TAMs.  The  role  of  TAMs  in  breast  tumor  tissue  is  not 
completely  understood.  In  the  tumor  microenvironment  TAMs  are  modified  to  secrete 
growth  factors  such  as  EGF,  VEGF,  TNFa  and  bFGF.^  It  has  been  suggested  that  TAMs 
play  a  role  in  neovascularization,  are  related  to  angiogenesis,  and  ultimately  survival  and 
proliferation  of  the  tumor  cell  population.® 

An  imaging  method  that  can  define  the  extent  of  this  activity  would  have  a 
significant  impact  on  cancer  imaging  and  research.  Since  macrophages  ingest  iron  oxide 
particles  we  investigated  the  application  of  optimized  imaging  methods  from  Phase  II  for 
measuring  iron  oxide  uptake  into  phagocytic  cells. 

In  vitro  experiments  have  demonstrated  our  ability  to  measure  iron  oxide  uptake 
into  a  macrophage-like  cell  line  and  to  characterize  uptake  efficiency  as  a  function  of 
particle  size,  dose  and  incubation  time.  The  results  of  these  experiments  indicate  that  (i) 
R2*is  more  than  100  times  more  sensitive  than  R2  to  iron  based  contrast  agents 
compartmentalized  in  cells,  (ii)  iron  uptake  by  phagocytic  cells  is  linear,  we  have  defined 
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a  parameter  for  uptake  efficiency  (dose*time)  which  allows  the  uptake  efficiency  to  be 
defined  for  any  cell  type  and  contrast  agent  (see  Figure  1),  and  (iii)  iron  within  cells  is 
quantifiable  using  R2'  mapping.  R2'  is  directly  proportional  to  iron  concentration  when 
compartmentalized  within  cells  (see  Figure  2).  This  data  is  included  in  the  submitted 
manuscript. 

In  Vivo  Animal  Studies  (Phase  IV) 

In  the  fourth  phase  of  this  research  project  MR  techniques,  developed  in  Phase  II 
and  testing  in  vitro  in  Phase  III,  are  applied  to  study  the  mammary  tumor  in  live  rats, 
theoretical  estimates  of  R2,  developed  in  Phase  II,  are  extended  to  predict  the  in  vivo 
detection  threshold  of  iron-oxide  in  cells. 

Images  were  acquired  of  both  mammary  tumor  tissue  and  normal  brain,  as  examples 
of  highly  homogeneous  and  heterogeneous  tissue  types,  respectively.  For  imaging  studies 
rats  were  anesthetized  with  an  intraperitoneal  injection  of  ketamine/xylazine  (100/10 
mg/kg)  and  positioned  supine  within  a  5  cm  diameter  solenoid  RF  coil  in  a  4  T 
Varian/Siemens  whole  body  scanner.  Axial  maps  of  both  rat  brain  and  tumor  were 
generated  using  the  3D  multi-echo  gradient  echo  sequence  (256x256, 32  slices  at  250  mm 
resolution,  6  cm  FOV,  77  kHz  full  receive  bandwidth,  120  ms  TR,  30°  flip  angle). 
Separate  image  volumes  were  produced  for  each  of  14  echoes  acquired  using  a  4  ms  echo 
spacing.  maps  were  acquired  using  a  single  slice  2D  multi-echo  CPMG  sequence 
(192x192,  4  mm  slice,  6  cm  FOV,  74  kHz  full  receive  bandwidth,  3  second  TR,  16 
echoes  with  7.5  ms  echo  spacing).  Separate  images,  zero  padded  to  256x256,  were 
produced  for  each  echo,  and  the  maximum  echo  time  was  acquired  at  TE,^  =  120  ms  (9.5 
minutes  per  acquisition).  Generation  of  R^*  and  R^  maps  from  the  corresponding  multi- 
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echo  gradient  echo  and  CPMG  images  was  accomplished  using  a  2  parameter  mono¬ 
exponential  fit  for  each  pixel.  The  single  slice  for  the  2D  CPMG  acquisition  was 
positioned  at  the  axial  center  of  the  3D  gradient  echo  volume  for  comparison  with  the 
maps.  R2  maps  were  generated  by  subtracting  the  R2  value  for  each  voxel  in  the  3D 
volume  from  the  R2  value  of  the  corresponding  in-plane  pixel  of  the  single  slice  “lap. 
While  this  method  has  some  disadvantages  due  to  partial  volume  effects  relative  to  the 
acquisition  of  a  full  3D  CPMG  volume,  the  approach  was  considered  adequate  since  little 
variation  in  the  axial  direction  was  observed  for  the  gradient  echo  images  along  the  4  mm 
range  of  the  2D  slice.  The  results  of  these  studies  are  shown  in  Figures  3-6. 
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Contrast  Agent  Development  (Phase  V) 

We  have  investigated  the  use  of  iron  oxide  nanoparticles  as  a  contrast  agent  for 
characterization  and  quantification  of  angiogenesis  in  mammary  tumors.  In  attempt  to 
further  increase  the  specificity  of  MR  contrast  agents,  superparamagnetic  particles  can  be 
linked  to  monoclonal  antibodies  (MAbs),  proteins  or  other  target-specific  moieties.  The 
original  goal  of  this  work  is  to  next  examine  iron  oxides  that  have  been  conjugated  to 
monoclonal  antibodies,  specifically  antibodies  directed  against  the  VEGF  receptor  that  is 
up  regulated  in  tumor  tissue.  We  have  successfully  developed  the  advanced  MRI 
techniques,  which  are  necessary  to  accurately  resolve  and  quantify  iron  content,  and  so 
are  well  suited  to  explore  the  next  phase  of  this  research.  We  have  developed  research 
collaboration  with  Dr.  Lee  Josephson,  a  world  expert  in  iron  oxide  based  contrast  agent 
development  at  the  Center  for  Molecular  imaging  research  at  MGH. 
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Key  Research  Accomplishments 

•  Implementation  and  characterization  of  a  carcinogen-induced  (ethyl-N-nitrosourea) 
breast  tumor  model  in  rats  that  produces  mammary  tumors  with  a  range  of 
malignancies,  from  benign  to  malignant. 

•  Implementation  and  optimization  of  MRI  pulse  sequence  technology  (multi-echo  3D 
gradient  echo  and  spin  echo  sequences)  for  iron  assessment. 

•  Design  and  coding  of  computer  programs  that  allow  automated  analysis  of  the  MR 
images  for  the  calculation  of  susceptibility  and  R2'  maps  for  determining  iron 
content. 

•  In  vitro  cell  studies:  Analysis  of  the  effect  of  particle  size,  incubation  time  and 
concentration  of  iron  oxide  nanoparticles  on  cellular  uptake  by  phagocytic  cells. 

•  In  vivo  animal  studies:  Analysis  of  the  particle  size,  dose  and  time  from  injection  of 
iron  oxide  nanoparticles  on  iron  uptake  in  TAMs. 

•  Correlative  analysis  of  MR  and  histologic  images 

The  stated  objective  of  developing  more  specific  contrast  agents  for  visualization 
of  tumor  microvasculature  has  not  been  accomplished  at  this  time.  This  phase  of  the 
research  has  not  begun  yet;  this  is  a  reflection  of  the  short  term  of  this  award.  We  are 
well  suited  to  move  forward  with  this  component  of  the  research,  however,  and  plan 
to  submit  a  grant  application  to  the  U.S.  Army  for  funds  to  continue  this  research. 
Intermediary  support  has  been  provided  by  The  Canadian  Breast  Cancer  Research 
Initiative  (see  Funding  below). 
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Reportable  Outcomes 

Presentations  and  Publications 

1 .  The  implementation  and  optimization  of  pulse  sequence  technology  for  iron 
assessment  was  presented  as  a  poster  at  the  9*  Meeting  of  the  International  Society 
for  Magnetic  Resonance  in  Medicine  (ISMRM:  Glasgow,  2001)  and  a  manuscript 
relating  to  this  work  has  been  accepted  for  publication  in  Magnetic  Resonance  in 
Medicine.  These  publications  are  appended. 

2.  A  second  manuscript  describing  the  application  of  these  methods  to  the  ENU  tumor 
model  is  in  preparation.  This  will  also  be  prepared  as  an  abstract  for  presentation  at 
the  10*  Meeting  of  the  International  Society  for  Magnetic  Resonance  in  Medicine 
(ISMRM;  Honolulu,  2002). 

Funding 

An  IDEA  Award  (one  year,  $50,000)  from  the  NCIC  (Canadian  Breast  Cancer 

Research  Initiative)  has  been  awarded  to  Dr.  Paula  Gareau  to  continue  the  research 

supported  by  this  grant. 
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Conclusions 

The  aim  of  this  Concept  Award  was  to  study  tumor  angiogenesis  in  breast  cancer, 
and  the  relationship  to  VEGF  and  TAMs,  using  advanced  MRI  techniques.  To  study 
these  features  of  breast  cancer  pathology  requires  a  method  that  is  noninvasive  and  will 
allow  frequent  and  repetitive  observations  in  the  same  tumor  over  time  during  tumor 
development.  Furthermore,  the  imaging  method  should  be  specific  (to  the  TAM  or 
VEGF  receptors),  sensitive  to  allow  accurate,  quantitative  measurements,  and  of  high- 
resolution  to  permit  visualization  of  the  microanatomical  structure  of  the  developing 
tumor. 

With  the  support  of  this  Concept  Award,  MRI  pulse  sequence  technology  has 
been  developed  and  optimized  for  the  relative  assessment  of  tissue  iron  concentrations 
based  on  the  quantitative  evaluation  of  T2  and  T2*  weighted  images.  These  methods 
allow  high-resolution  MR  images  to  be  obtained  in  a  3D  format.  Computer  programs 
were  written  for  the  automated  analysis  of  the  MR  images  for  the  calculation  of 
susceptibility  and  R2'  maps  for  determining  iron  content.  These  pulse  sequence  and 
software  developments  were  validated  in  vitro  in  a  cell  preparations  and  in  vivo  in  an 
ENU-induced  mammary  tumor  model.  We  have  defined  a  parameter  for  uptake 
efficiency  (dose*time)  that  allows  the  uptake  efficiency  to  be  defined  for  any  cell  type 
and  contrast  agent.  Iron  oxide  nanoparticles  were  studied  for  use  as  MRI  contrast  agents 
for  imaging  TAMs  in  mammary  tissue.  We  have  successfully  developed  MRI  techniques, 
which  are  necessary  to  accurately  resolve  and  quantify  iron  content  in  vivo  and  are 
prepared  to  investigate  the  next  phase  of  this  research. 
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Assessing  the  role  of  the  TAM  in  tumor  angiogenesis  and  tumorigenesis  using 
this  approach  may  result  in  novel  in  vivo  markers  for  tumor  type,  aggressiveness  and 
tumor  margin  delineation.  Furthermore,  the  association  of  novel  contrast  agents  with 
VEGF  receptors,  TAMs  and/or  phagocytosis  of  contrast  agents  by  TAMs,  forms  the  basis 
of  a  potential  therapeutic  strategy  in  which  the  contrast  agent  could  serve  as  a  targeted 
carrier  of  drugs. 
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Figure  1:  Cellular  uptake  of  two  different  sized  iron  oxide  contrast  agents  iResovist® 
(60nm)  and  Resovist-S®  (20nm).  The  slope  of  the  line  represents  the  uptake 
efficiency  (UE).  Uptake  by  phagocytosis  is  linear  so  a  unique  uptake  efficiency  is 
defined  for  every  cell  and  contrast  agent 
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Figure  2;  Rj'  relaxivity  at  4  Tesla.  Experiments  with  different  sized  iron  oxide 
contrast  agents  (Resovist=60nm,  Resovist-S=20nm)  in  phagocytic  cells  or  free  in 
solution.  R2'  is  independent  of  contrast  agent  (or  field  strength)  when 
compartmentalized  in  cells. 


Figure  3:  Example  gradient  echo  images  and  R2  maps  from  rat  brain  (a,b)  and  tumor  (c,d).  3D 
multi-echo  gradient  echo  images  are  of  the  first  echo  (4.5  ms  TE).  Dashed  square  indicates 
typical  ROI  selected  for  background  R2'  analysis.  Grayscale  bars  (0  to  200  s'^)  for  R2'  maps 
(b,d)  are  identical. 
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Figure  4:  Summary  of  in  vivo  background  distribution  parameters  for  rat  tumor  and 
brain.  Bar  chart  indicates  background  R2  mean  {Rziyg),  standard  deviation  (a^p,  median, 
and  the  threshoid  greater  then  95%  {pgg)  and  99%  (p,,,)  of  all  values  within  the 
ROI.  Data  is  averaged  for  3  animals  and  error  bars  indicate  the  95%  confidence  interval. 
Small  arrows  indicate  normal  distribution  predictions  of  Pgg  and  p^,  based  on  the 
reported  mean  {Rs^g)  and  standard  deviation  (0^^).  Agreement  between  the  mean  and 
median  values  reflects  distribution  symmetry.  Agreement  between  the  arrows  and  the 
measured  p  values  indicates  normally  distributed  Revalues. 
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Figure  5:  Method  for  separating  mesoscopic  {R2mes')  macroscopic  {R2mac)  components  of 
background  R2  {R2bg')-  Maps  of  R2mes'  (3)  and  R2mac'  (b)  from  rat  tumor  generated  using  the 
method  illustrated  in  (c).  The  plot  (c)  shows  relaxation  rate  dependence  on  pixel  size  for  a 
single  pixel  from  a  portion  of  the  tumor  having  a  large  macroscopic  R2'  component  (white 
arrows  in  (a)  and  (b)).  The  resolution  dependence  of  R2  (circles)  and  R2  (squares)  is  shown 
along  with  the  linear  fit  for  each  (dashed  and  dotted  line,  respectively).  R2bg'  (solid  line)  is 
derived  from  subtraction  of  the  R2  and  R2  fitted  curves.  R2mes  ^f^^l  R2mac  the  gradient  echo 
acquisition  resolution  (Ax^,)  is  indicated  (see  Eq.  [3.3]).  Maps  of  R2mac  (h)  have  little  fine 
structure  and  correlate  with  regions  of  large  scale  field  inhomogeneity  while  R2mes'  maps  (a) 
correlate  with  local  tissue  properties. 
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Figure  5:  Summary  of  rat  brain  and  tumor  relaxation  rates  separated  into  mesoscopic 
and  macroscopic  components.  The  mesoscopic  (dotted  bar)  and  macroscopic  (angled 
line  bar)  components  of  R2’,  R2  and  f?2bg'.  measured  using  the  method  described  in  Rg.  4, 
are  shown  for  the  brain  (a)  and  tumor  (b)  ROI.  Note  the  different  scale  on  the  axis.  Data  is 
averaged  for  3  animals  and  error  bars  indicate  the  95%  confidence  interval.  The  fractional 
component  of  background  Rsif^ac)  gradient  echo  acquisition  resolution  (AXo=  234 
pm)  is  33%  and  25%  for  rat  brain  and  tumor,  respectively. 
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Abstract 


The  relaxation  rates  of  iron-oxide  nanoparticles  compartmentalized  within  cells  were  studied 
and  found  to  satisfy  predictions  of  static  dephasing  (SD)  regime  theory.  THP-1  cells  in  cell 
culture  were  loaded  using  2  different  iron-oxide  nanoparticles  (SPIO  and  USPIO)  with  4 
different  iron  concentrations  (0.05,  0.1,  0.2,  0.3  mg/ml)  and  for  5  different  incubation  times  (6, 
12,  24,  36,  38  hours).  Cellular  iron-oxide  uptake  was  assessed  using  a  newly  developed  imaging 
version  of  MR  susceptometry  and  found  to  be  linear  with  both  dose  and  incubation  time. 
sensitivity  to  iron-oxide  loaded  cells  was  found  to  be  70  times  greater  than  for  R2  and  3100  times 
greater  than  for  Ry.  This  differs  greatly  from  uniformly  distributed  nanoparticles  and  is  consistent 
with  a  cellular  bulk  magnetic  susceptibility  relaxation  mechanism.  Cellular  magnetic  moment 
was  large  enough  that  72;'  relaxivity  agreed  closely  with  SD  regime  theory  predictions  for  all  cell 
samples  tested.  Uniform  suspensions  of  SPIO  and  USPIO  produced  R2  relaxivities  that  were  a 
factor  of  3  and  8  less  than  SD  regime  theory  predictions,  respectively.  These  results  are 
consistent  with  theoretical  estimates  of  the  required  mass  of  iron  per  compartment  needed  to 
guarantee  SD-regime-dominant  relaxivity.  For  cellular  samples,  R2  was  shown  to  be  dependent 
on  both  the  concentration  and  distribution  of  iron-oxide  particles,  while  was  sensitive  to  iron- 
oxide  concentration  alone,  being  independent  of  both  field  strength  and  particle  size.  This  work 
is  an  important  first  step  in  quantifying  cellular  iron  content  and  ultimately  mapping  the  density 
of  a  targeted  cell  population. 

Keywords: 

Inflammation,  Susceptibility,  Contrast  Agent,  Iron-Oxide,  Relaxometry 
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Introduction 


The  macrophage  is  a  white  blood  cell  of  the  immune  system  which  is  intimately  involved  in 
the  development  and  progression  of  inflammatory  disease  (1).  Knowledge  of  the  time  course, 
spatial  distribution  and  magnitude  of  inflammatory  response,  as  well  as  the  evolution  of 
macrophage  activation  state,  is  critical  for  understanding  disease  progression  (2,3).  Currently,  the 
most  common  method  for  assessing  the  extent  of  inflammation,  or  the  effectiveness  of  various 
therapeutic  strategies  that  inhibit  inflammation,  involves  a  histopathologic  evaluation  of  multiple 
tissue  sections.  Imaging  methods  that  are  capable  of  defining  the  extent  of  inflammatory  activity 
in  vivo  would  have  a  significant  impact  on  inflammation  research. 

It  has  recently  become  possible  to  visualize  macrophage  cells  in  vivo.  Weissleder  and 
colleagues  have  pioneered  this  development,  utilizing  iron  oxide  nanoparticles  to  image  cellular 
uptake  and  trafficking  by  MRI  (4).  This  technique  exploits  the  ability  of  certain  cell  types  to 
ingest  small  particles  in  culture  and  in  vivo  through  phagocytosis.  For  example,  intravenously 
administered  iron  oxide  particles  accumulate  in  macrophages  of  acute  lesions  in  experimental 
autoimmune  encephalomyelitis  (5),  and  in  tumor  cells  (6).  While  many  papers  report  sensitivity 
to  the  presence  of  iron-oxide  using  T,  (6,7),  T2  (5,8)  and  T2  (9,10)  weighted  imaging,  methods 
quantifying  the  concentration  of  iron-oxide  in  vivo,  as  a  first  step  towards  characterizing  the 
degree  of  inflammatory  cell  infiltration,  are  still  lacking.  The  question  of  which  acquisition 
strategy  is  the  most  sensitive  and,  in  particular,  the  most  accurate  for  quantifying  the  amount  of 
iron-oxide  compartmentalized  within  cells  is,  thus,  still  an  open  one.  We  have  investigated,  both 
theoretically  and  experimentally,  the  relaxation  mechanism  responsible  for  NMR  signal  decay  in 
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the  presence  of  iron-oxide  nanoparticles,  both  in  uniform  suspension  and  following  uptake  into 
cells.  Our  results  represent  an  important  first  step  toward  quantifying  the  cellular  iron  content 
and  ultimately  mapping  the  density  of  the  targeted  cell  population. 

Theory 

Iron-oxide  nanoparticles  are  superparamagnetic  (SPM)  particles  (11),  with  magnetic 
moments  more  than  3  orders  of  magnitude  greater  than  those  of  paramagnetic  contrast  agents 
(12).  The  quantitative  dependence  of  both  Rj  and  R2  (inverse  Tj  and  T^)  has  been  estimated  using 
the  Solomon-Bloembergen-Morgan  equations,  modified  for  the  case  of  SPM  particles,  by 
Koenig  (13),  Bulte  (12)  and  others.  Below  1  T,  these  traditional  quantum  solutions  describing 
electron-proton  dipolar  relaxation  are  quite  successful  in  quantifying  Rj  and  R2  relaxation  rates 
for  iron-oxide  particles  suspended  in  gel.  SPM  particles  are  observed  to  have  a  much  smaller 
effect  on  Tj  then  on  relaxation  compared  to  the  effects  produced  by  paramagnetic  agents, 
owing  to  the  large  SPM  magnetic  moment. 

However,  when  SPM  particles  are  compartmentalized  within  cells  rather  than  uniformly 
suspended  in  gel,  the  quantum  solution  fails.  This  effect  has  been  observed  both  in  vivo,  for  SPM 
particles  sequestered  within  Kupffer  cells  of  the  reticuloendothelial  system  of  the  liver  (14),  and 
in  vitro,  in  cell  culture  (4).  In  these  studies,  the  ratio  of  R2'  to  R2  is  seen  to  increase  dramatically 
compared  to  that  observed  in  gel  suspension,  where  the  quantum  solution  would  predict  identical 
values  for  R2*  and  R  2.  The  quantum  solution  relies  on  the  extreme  motional  narrowing 
approximation,  which  states  that  water  diffusion  between  SPM  particles  occurs  on  a  much 
shorter  time  scale  than  the  peak  resonant  frequency  offset.  When  iron-oxide  particles  are 
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compartmentalized,  as  within  cells  for  instance,  the  cellular  scale  bulk  magnetic  susceptibility 
(BMS)  effects  produce  dipolar  field  inhomogeneities  large  enough  to  violate  this  assumption. 
Unfortunately,  no  analytic  solution  describing  Rj  or  R2  exists  in  this  large  perturber  case  (15), 
though  Monte  Carlo  simulations  of  water  diffusing  through  the  BMS  induced  dipolar  fields  are 
quite  successful  in  predicting  R2  and  R2'  (16-20). 

Monte  Carlo  predictions,  though  useful  for  understanding  the  dominant  relaxation 
mechanism,  are  difficult  to  generalize  for  differences  in  compartment  size,  field  strength  and 
pulse  sequence  parameters  such  as  refocusing  pulse  spacing.  However,  an  analytic  solution  does 
exist  when  perturber  strength  becomes  large  enough  to  satisfy  the  static  dephasing  (SD)  regime 
condition  (21,22).  In  this  regime,  the  large  compartmental  magnetic  moment  produces  a  strong 
enough  outer  sphere  dipolar  field  such  that  diffusion  has  no  effect  on  the  NMR  signal  decay.  In 
the  limit  of  the  SD  regime,  Rj  and  R2  enhancement  due  to  the  presence  of  SPM  particles 
approaches  zero  (15),  while  R2',  as  predicted  by  Yablonskiy  and  Haacke  (22)  in  Eq.  [58]  of  their 
work,  is  given  by  the  following  relation  (converted  to  SI  units): 

where  y  is  the  proton  gyromagnetic  ratio  and  LMD  is  the  local  magnetic  dose,  defined  by 
Yablonskiy  and  Haacke,  as  the  tissue  averaged  magnetization  difference  due  to  the  presence  of 
SPM  particles.  By  this  formalism,  LMD  =  where X  is  the  fractional  volume  occupied  by  Ihe 

magnetic  compartments  (cell  volume  fraction  for  instance)  and  is  the  internal  magnetization 
of  the  compartment  due  to  the  presence  of  SPM  particles.  R2'  is  defined  as  R2*  -  R2.  In  this  case. 
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R2  reflects  the  enhanced  relaxation  rate  observed  with  gradient  echo  acquisitions  due  to  the 
presence  of  SPM  particles. 


The  conditions  for  which  the  outer  sphere  static  dephasing  regime  equation,  Eq.  [1], 
accurately  predicts  R2  relaxation  is  also  indicated  by  Yablonskiy  and  Haacke.  This  SD  regime 
limit  is  given  in  Eq.  [55]  of  their  work  (22)  as  follows  (converted  to  SI  units): 


R:  » 


6D 


8(0,  = 


[2] 


where  D  is  the  diffusion  coefficient,  R„  is  the  mean  magnetic  compartment  size  and  5t0s  is  the 
resonant  frequency  offset  at  the  external  edge  on  the  equator  of  the  magnetic  compartment, 
called  the  equatorial  frequency  offset.  This  SD  regime  limit  definition  is  somewhat  non- 
traditional  as  it  requires  a  large  diffusion  time  from  one  magnetic  compartment  location  to 
another  relative  to  the  inverse  frequency  offset  midway  between  particles.  This  differs  from  the 
traditional,  and  less  restrictive,  definition  requiring  diffusion  time  across  the  particle  dimension, 
R„  to  be  large  relative  to  the  inverse  equatorial  frequency  offset,  bcOs"'. 

Though  well  described  by  Yablonskiy  and  Haacke,  the  fundamental  theory  provided  by  Eqs. 
[1]  and  [2]  has  not  been  confirmed  experimentally  to  date.  The  authors  suggest  application  of 
their  theory  to  quantifying  free,  uncompartmentalized  iron-oxide  nanoparticles,  however,  the 
static  dephasing  regime  equations  are  likely  to  be  more  suited  to  the  application  of  iron-oxide 
compartmentalized  within  cells.  In  this  work,  the  validity  of  predictions  for  R2’  as  well  as  the  SD 
regime  limit  are  tested  for  both  compartmentalized  and  free  iron-oxide  suspensions. 


6 


Methods 


THP-l  Cell  Preparation 

THP-1  cells  were  chosen  to  study  intracellular  iron  signal  effects  as  they  can  be  prepared 
with  a  wide  range  of  intracellular  iron  concentrations  through  phagocytosis  of  iron  core  SPM 
contrast  agents.  THP-1  cells  are  a  human  promonocytic  cell  line  derived  from  acute  monocjdic 
leukemia  (from  the  American  Type  Culture  Collection).  The  cells  were  grown  in  25-ml  flasks 
with  RPMI-1640  medium  (Gibco)  supplemented  with  10,000  units/ml  penicillin,  10  mg/ml 
streptomycin  and  10%  FBS  at  37°C.  5  x  10®  cells  placed  in  5  ml  medium/dish  were  activated  by 
adding  100  ng/ml  of  Phorbol  12-myristate  13-acetate  for  5  days.  After  activation,  non-adhered 
cells  were  removed  through  medium  changes  leaving  approximately  4x10®  cells/dish. 

Two  carboxydextran-coated  magnetite  nanoparticles  were  compared:  SHU  555A  (Schering 
AG,  Berlin),  a  superparamagnetic  iron  oxide  (SPIO)  contrast  agent  with  a  mean  particle  size  of 
about  60  nm;  and  a  new  ultrasmall  superparamagnetic  iron  oxide  (USPIO)  contrast  agent,  SHU 
555C  (Schering  AG,  Berlin),  having  a  mean  particle  size  of  about  20  nm.  Iron  loading  was 
accomplished  by  adding  either  SPIO  or  USPIO  contrast  agent  at  one  of  4  doses  (0.05,  0.1,  0.2, 
0.3  mgFe/ml)  and  5  incubation  times  (6, 12, 24,  36, 48  hours).  The  dose  group  was  incubated  for 
24  hours  and  the  incubation  time  group  using  a  dose  of  0.1  mgFe/ml.  After  iron  incubation,  the 
adhered  cells  were  detached  by  trypsinization  (adding  0.05%  trypsin  for  3  minutes  on  ice)  and 
washed  in  HBSS  three  times.  The  cell  viability  was  then  confirmed  by  the  trypan-blue  exclusion 
test  and  the  number  of  cells  were  counted. 

NMR  samples  for  the  dose  and  incubation  time  group  were  prepared  with  a  cellular 
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concentration  of  4  x  10®  cells/ml.  This  corresponds  to  a  cellular  volume  fraction  of  only  0.189%, 
given  the  reported  THP-1  cell  volume  (23)  of  473  ±  16  |im^.  Each  1  ml  sample  was  placed  within 
5  mm  NMR  tubes  and  suspended  in  2%  w/w  gelatin  to  prevent  sedimentation.  A  cell  density 
group,  incubated  for  24  hours  at  0.1  mgFe/ml,  was  prepared  by  concentrating  cells  from  multiple 
culture  dishes  after  iron-oxide  incubation.  1  ml  NMR  samples  were  prepared  at  5  cellular 
densities  (1,  2,  4,  8,  16  x  10®  cells/ml).  In  total,  4  dose  group  samples,  4  unique  incubation  time 
samples,  5  cell  density  samples  and  a  control  sample,  incubated  for  24  hours  without  contrast 
agent  present,  were  prepared  using  each  contrast  agent  (see  table  1). 

Cellular  Iron  Quantification 

Two  properties  of  cellular  iron  loading  were  monitored:  intracellular  magnetization  and  iron 
mass  per  cell.  Direct  measurement  of  intracellular  magnetization  is  important  since  this  is  the 
property  that  appears  in  theoretical  treatments  of  susceptibility-induced  NMR  signal  response. 
Magnetization  was  measured  using  a  susceptometry  technique,  the  Reilly-McConnell- 
Meisenheimer  (RMM)  method  (24,25),  modified  in  this  work  to  use  multi-echo  imaging. 
Compared  to  measurements  of  spectral  broadening,  as  employed  by  the  original  RMM  method, 
this  approach  more  accurately  characterizes  the  susceptibility-induced  field  inhomogeneity 
pattern  outside  of  a  NMR  tube  oriented  perpendicular  to  B„.  This  NMR  based  susceptometry 
method  permits  non-destructive  measurement  of  sample  magnetization  even  in  situations  where 
extremely  short  relaxation  times  preclude  conventional  NMR  relaxometry  or  spectroscopy. 

Bulk  magnetic  susceptibility  for  any  substance  within  a  NMR  tube  can  be  estimated  from  the 
magnetic  field  inhomogeneity  pattern  established  outside  such  a  NMR  tube  oriented  orthogonal 
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to  B„.  The  solution  to  Maxwell’s  equation  for  the  field  offset  produced  external  to  co-axial 
cylindrical  boundaries  of  differing  susceptibility  (glass  walled  NMR  tube)  is  (24): 


AB(r)  = 


B„  cos 20 


3 


]  [^{X2-Xi)  +  alix^-Xi)] 


[3] 


where  Xa  atid  X3,  are  the  bulk  magnetic  susceptibility  within  the  NMR  tube  (desired),  of  the 
NMR  glass  and  of  the  material  external  to  the  NMR  tube  (distilled  water),  respectively,  a,  and 
are  the  internal  and  external  radii  of  the  NMR  tube  (2.1,  2.5  mm),  respectively,  r  is  the  radius  to 
the  field  profile  position  of  interest  external  to  the  NMR  tube  (i.e.  r  =  |r|  where  r  is  the  vector 
from  tube  origin  to  position  of  interest)  and  0  is  the  angle  between  r  and  B„.  For  samples  with 
saturated  magnetization,  magnetization  does  not  scale  with  static  field  and  x,5o  is  replaced  by 
LMD,  the  local  magnetic  dose  (22).  This  is  the  case  at  normal  MRI  field  strengths  for 
superparamagnetic  species  like  the  two  iron  oxide  nanoparticles  used  in  this  work  (1 1,12,22). 

The  RMM  apparatus  consisted  of  a  5  mm  outer  diameter  NMR  tube  concentrically  placed, 
using  teflon  spacers,  within  a  15  mm  inner  diameter,  10-ml  graduated  cylinder.  The  graduated 
cylinder  contained  distilled  water  and  the  surfactant,  Photo-Flo  200  (0.5%),  to  inhibit  bubble 
formation.  The  apparatus  was  placed  vertically  within  a  5  cm  diameter  solenoid  RF  coil  in  a  4  T 
Varian  whole  body  scanner.  The  field  inhomogeneity  pattern  was  mapped  by  generating  a  field 
map  from  the  phase  evolution  versus  echo  time  in  images  acquired  using  a  2D  multi-echo 
gradient  echo  sequence  (128  x  128, 10  mm  slice,  6  cm  FOV,  8  ms  echo  spacing,  20  echoes).  The 
field  map  was  computed  through  a  Fourier  transform  of  the  complex  signal  along  the  echo  train 
for  each  pixel.  The  peak  frequency  for  each  pixel  from  the  magnitude  of  the  Fourier  transform 
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defined  the  local  field  offset.  A  sample  field  map  (0.05  mg/ml  SPIO  in  gel)  is  shown  in  Fig.  lb. 
The  corresponding  magnitude  image  is  shown  in  Fig.  la.  Phase  wrap  due  to  aliasing  occurs  for 
field  offsets  greater  than  the  Nyquist  frequency  {+\IQ.»esp)  or  ±62.5  Hz). 

Field  map  inhomogeneity  external  to  the  NMR  tube  was  characterized  by  averaging  4  field 
profiles  extending  radially  outward  from  the  external  edge  of  the  NMR  tube  in  each  orthogonal 
direction.  Figure  Ic  shows  the  4  profiles  corresponding  to  the  field  map  of  Fig.  lb.  To  correct  for 
aliasing,  each  profile  was  imwrapped  along  the  radial  direction  for  frequency  discontinuities 
greater  than  half  the  Nyquist  frequency.  Useful  field  profiles  could  thus  be  obtained  as  long  as 
the  difference  in  frequency  offsets  of  adjacent  pixels  did  not  exceed  the  Nyquist  frequency. 
Unwrapping  discontinuities  associated  with  aliasing  for  an  entire  image  is  challenging,  so  the 
inhomogeneity  pattern  was  characterized  using  only  4  profiles  to  avoid  these  difficulties.  To 
enhance  SNR  of  the  profiles  and  counter  the  effect  of  any  global  linear  field  variation  across  the 
image,  the  4  profiles  were  averaged  together  (after  appropriate  polarity  reversal  of  the  two  0  = 
90°  profiles). 

This  average  profile  was  fit  to  the  theoretical  inhomogeneity  pattern,  shown  in  Eq.  [3],  to 
determine  the  susceptibility  within  the  NMR  tube.  A  two  parameter  constrained  least  squares  fit 
of  Eq.  [3]  was  used.  The  average  profile,  along  with  the  results  of  the  fit,  are  shown  in  Fig.  Id. 
The  two  free  parameters  in  the  fit  were  Xj,  the  sample  susceptibility,  and  a  baseline  constant 
added  to  account  for  global  resonance  offset.  The  susceptibility  of  water  (Xj)  and  NMR  tube 
glass  (X2)  were  assumed  to  be  -9.06  ppm  and  -9.77  ppm,  respectively.  X2  was  measured  using  a 
water-filled  tube,  and  fitting  Eq.  [3]  for  X2.  with  Xi  and  Xs  constrained  to  equal  the  susceptibility 
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of  water.  Sample  susceptibility  was  converted  to  mg/ml  of  Fe  through  calibration  with  known 
concentration  reference  samples  of  SPIO  and  USPIO,  in  the  form  of  uniform  suspensions  of  the 
nanoparticles  in  gelatin. 

Iron  concentrations  for  several  samples  were  measured  using  atomic  emission  spectroscopy 
(AES).  This  method  directly  quantifies  iron  mass  independent  of  valence.  AES  is  also  readily 
available  in  most  hospital  labs  and  hence  represents  a  practical  gold  standard  for  validating  the 
susceptometry  method  estimate  of  iron  concentration.  AES  measurements  were  also  used  to 
verify  that  no  cell  lysis  in  cellular  samples  had  occurred  since  lysis  would  cause  the  undesired 
release  of  iron-oxide  into  the  extra-cellular  gelatin.  This  was  accomplished  by  gently  heating  the 
samples  to  37°C  (after  all  NMR  measurements),  separating  the  gelatin  from  the  cells  via 
centrifuging  and  performing  AES  on  the  isolated  gelatin.  On  average  for  the  3  samples  tested, 
only  6%  of  the  iron  in  cellular  samples  was  extracellular;  this  experiment  confirmed  that 
extracellular  iron  could  be  ignored  in  our  analysis. 

Relaxation  Mechanism 

The  methods  described  in  this  section  serve  three  purposes:  to  identify  the  dominant 
relaxation  mechanism  affecting  NMR  signal  for  iron-oxide  loaded  cells;  to  measure  the  relative 
sensitivity  of  different  relaxation  rates  {R„  R2  and  and  to  test  the  accuracy  of  static 
dephasing  regime  predictions  in  this  case.  Measurements  were  made  at  1.89  T  on  a  SMIS/Oxford 
30  cm  bore  system  and  at  4  T  on  a  Varian/Siemens  whole  body  scanner. 

Measurements  were  made  at  1.89  T  using  a  2D  hybrid  relaxometry  pulse  sequence  which 
simultaneously  provided  T,  and  T2  relaxometry  spectra  (26).  This  approach  produces  both  T;  and 
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results  in  a  scan  time  of  approximately  5  minutes  per  sample  (described  in  greater  detail  in 
Ref.  26).  The  experimental  arrangement  consisted  of  sample  NMR  tubes  being  placed  within  a 
15  mm  diameter  solenoid  RF  coil  orthogonal  to  the  main  field.  Acquisition  consisted  of  a  non- 
selective  CPMG  train  using  a  90-180  pulse  spacing,  T^p^g,  of  1.25  ms.  The  decay  curve  was 
sampled  with  950  even  echo  times  spanning  5  to  4750  ms.  An  inversion  recovery  pulse  preceded 
the  CPMG  train  with  30  different  inversion  times,  77,  spanning  40  to  7000  ms.  Processing 
consisted  of  a  2D  non-negative  least  squares  (NNLS)  fit  using  200  bins  spanning  0-2000  ms 
and  35  T,  bins  spanning  0-3500  ms.  The  2D  spectra  for  all  samples  were  dominated  by  a  single 
peak  and  so  the  spectra  were  collapsed  and  reported  as  a  single  value  of  Tj  and  corresponding 
to  that  peak.  T/  was  estimated  using  this  same  experimental  arrangement  at  both  1.89  T  and  4  T 
through  measurement  of  the  spectral  linewidth  from  a  pulse-acquire  sequence.  The  7?^’  relaxation 
rate  was  estimated  as  7?^*  =  n-FWHM  where  FWHM  is  the  full  width  at  half  maximum  of  the 
Lorentzian  spectral  line  in  Hz. 

In  order  to  further  study  the  relaxation  mechanism  and  test  the  static  dephasing  regime 
approximation,  7^  measurements  were  made  at  4  T  using  a  2D  multi-echo  CPMG  imaging  pulse 
sequence  (128  x  128,  10  mm  slice,  6  cm  FOV).  Other  than  the  pulse  sequence,  the  experimental 
arrangement  was  identical  to  that  described  for  NMR  susceptometry  in  the  previous  section. 
Images  were  reconstructed  for  every  echo  of  the  CPMG  train.  Two  different  90-180  CPMG 
pulse  intervals,  x^p^g,  were  used  in  separate  acquisitions  to  sample  200  ms  of  the  echo  train  (20 
echoes  with  x„p„g  =  5  ms  and  10  echoes  with  x^g  =  10  ms).  The  signal  decay  curve,  generated 
from  a  region  of  interest  within  the  NMR  tube,  was  fit  using  ID  NNLS  to  produce  a  7^  spectrum. 
Only  one  peak  was  found  for  all  samples  and  the  peak  7^  was  used  to  report  the  sample  7^  value. 
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To  ensure  that  the  estimated  was  not  biased  by  the  different  number  of  echo  times  sampling 
the  signal  decay  curve  for  each  X(,p„g,  the  =  5  ms  signal  decay  curve  >vas  separately  fit  using 
either  all  echoes  or  only  even  echoes.  No  significant  differences  were  observed. 

Results 

Cellular  Iron  Quantification 

The  susceptometry  calibration  curves  for  the  SPIO  and  USPIO  iron-oxide  contrast  agents  are 
shown  in  Fig.  2.  The  slope  of  the  calibration  curves  are  1319  ±  69  and  1360  ±  55  mG  per  mg/ml 
of  Fe  for  SPIO  and  USPIO,  respectively.  All  reported  uncertainties  refer  to  95%  confidence 
intervals  throughout  this  section.  The  two  SPM  particle  calibration  slopes  agree  within 
uncertainty  (p=0.05).  For  single  magnetic  domain,  superparamagnetic  iron-oxide  particles  having 
the  same  iron  spin  state,  this  should  be  the  case  (11,12).  Since  both  SPM  particles  studied  are 
preparations  of  similar  4.6  nM  magnetite  cores,  differing  only  in  number  of  cores  per  particle 
and  overall  particle  size  (27),  agreement  between  the  calibration  slopes  is  not  unexpected. 

Table  1  summarizes  the  cellular  iron  uptake  data  for  SPIO  and  USPIO.  The  results  for  the 
dose,  incubation  time  and  cellular  density  groups,  as  described  in  the  methods  section,  are 
shown.  AES  iron  concentration  for  3  SPIO  and  3  USPIO  samples  are  shown  in  curly  brackets  for 
comparison  with  the  susceptometry  method  and  agree  within  xmcertainty  (p=0.05)  with  a  7%  rms 
error.  Figure  3  plots  iron  uptake  in  pgFe/cell  for  SPIO  and  USPIO.  For  each  contrast  agent,  iron 
uptake  for  both  dose  and  incubation  time  groups  was  plotted  collectively  versus  the  incubation 
time  -  dose  product.  As  expected,  higher  doses  and  longer  incubation  times  resulted  in  increased 


13 


uptake.  This  uptake  was  linear  when  plotted  against  either  dose  or  incubation  time  independently 
(not  shown)  with  all  curves  having  correlation  coefficients  greater  than  0.97. 

Assuming  this  linearity  remains  valid  for  a  distribution  of  doses  and  incubation  times  aroxmd 
the  reference  dose  and  incubation  time  (0.1  mgFe/ml  and  24  hours),  an  uptake  efficiency  can  be 
defined  about  this  reference.  In  Fig.  3,  uptake  efficiency  is  reflected  by  the  slope  of  iron  uptake 
(pgFe/cell)  versus  the  product  of  incubation  time  and  dose,  using  data  from  both  dose  and 
incubation  time  groups.  The  correlation  coefficient  was  0.97  and  0.98  for  the  SPIO  and  USPIO 
uptake  curves,  respectively,  supporting  the  legitimacy  of  a  linear  fit.  The  uptake  efficiencies  for 
SPIO  and  USPIO  were  determined  to  be  3.06  ±  0.40  and  0.80  ±  0.08  pgFe/cell  per  hour»dose, 
respectively.  There  was  no  evidence  of  iron  uptake  saturation,  or  a  reduced  rate  of  iron  uptake 
after  time,  during  phagocytosis  for  THP-1  cells  over  the  range  of  doses  and  incubation  times 
studied.  This  is  perhaps  not  surprising  considering  that  phagoc)dosis  is  a  concentration  mediated 
process  that  eventually  results  in  cell  death  rather  than  saturation.  Our  cell  viability  tests 
indicated  that  even  at  the  highest  uptake,  cell  death  had  not  occurred  to  any  significant  extent. 

Relaxation  Mechanism 

Figure  4  shows  the  relative  sensitivity  of  different  NMR  relaxation  rates  at  1.89  T  to  the 
presence  of  iron-oxide  compartmentalized  within  THP-1  cells.  The  relaxation  rates  are  plotted 
for  the  constant  cellular  density  SPIO  and  USPIO  groups.  A  log-log  plot  was  used  because  of 
variation  over  2  orders  of  magnitude  in  the  R„  R^  and  R2'  relaxation  rates.  Table  2  lists  the 
relaxivity  of  each  NMR  time  constant  obtained  from  linear  regression  for  both  SPIO  and  USPIO 
in  cells.  In  each  of  the  7?;,  R2  and  7?^*  cases,  relaxivities  for  SPIO  and  USPIO  are  similar,  agreeing 
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to  within  40%,  although  different  at  the  95%  confidence  level.  On  average  for  each  iron-oxide 
particle,  R2  is  70  times  greater  than  R2  and  3100  times  greater  than  Ri,  for  SPM  loaded  cells.  As 
mentioned  in  the  Theory  section,  this  observation  is  consistent  with  cellular  scale  bulk  magnetic 
susceptibility  effects  violating  the  extreme  motional  narrowing  assumption  implicit  in  quantum 
dipolar  relaxation  theory  (13,15).  However,  the  question  of  whether  the  cellular  magnetic 
moment  is  large  enough  to  satisfy  conditions  of  SD  regime  theory,  described  by  Yablonskiy  and 
Haacke  (22),  still  remains. 

Figure  5  demonstrates  the  impact  of  iron-oxide  compartmentalization  within  cells  on  R2 
relaxivity.  The  R2  relaxation  rate  curves,  measured  at  4  T,  are  plotted  for  SPIO  and  USPIO  in 
gel  suspensions  as  well  as  compartmentalized  within  cells.  All  data  from  the  dose,  incubation 
time  and  cellular  density  groups  are  included.  It  is  important  to  recall  that  each  value  of  LMD 
corresponds  to  a  specific  iron-oxide  concentration  so  that  differences  between  curves  reflect 
differences  due  to  the  microscopic  distribution,  not  to  the  concentration  of  iron-oxide.  The  static 
dephasing  regime  theoretical  prediction  for  R2',  described  in  Eq.  [1],  is  also  plotted  for 
comparison.  The  intercept  for  this  curve  is  the  R2  relaxation  rate  measured  from  the  control  cell 
suspension  that  was  not  exposed  to  iron-oxide  nanoparticles.  Although  R2  is  given  by  the 
addition  of  R2  and  R2,  the  R2  component  was  neglected  since  table  2  indicates  that  R2  is  at  least 
70  times  less  than  ^nd  hence  R2.  The  theory  curve  represents  the  maximum  relaxation  rate 
possible  as  described  by  the  static  dephasing  regime  theory. 

The  R2*  relaxation  rates  per  unit  LMD  (“relaxivity”)  obtained  from  the  linear  regression 
slopes  of  each  curve  are  listed  in  table  3.  The  relaxivities  for  SPIO  at  4  T  and  USPIO  at  1.89  T 
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agree  with  SD  regime  theory  within  the  95%  confidence  interval,  while  the  relaxivities  for  SPIO 
at  1.89  T  and  USPIO  at  4  T  are  40%  higher.  Outer  sphere  dipolar  field  SD  regime  theory  predicts 
a  relaxivity  of  10.78  s"'  mG''  (Eq.  [1]).  No  field  strength  dependence  in  R2*  relaxivity  was 
observed:  these  relaxivities  averaged  across  the  two  iron-oxide  agents  compartmentalized  within 
cells  were  13.5  ±  2.3  and  12.8  ±  1.6  s"'  mG  ’  at  1.89  T  and  4  T,  respectively.  This  is  expected  for 
iron-oxide  contrast  agents  since  saturation  of  SPM  magnetization  occurs  for  any  field  strength 
greater  than  1  T  (1 1,12,22).  Similarly,  no  dependence  on  SPM  particle  size  was  observed:  the  R2* 
relaxivities  averaged  over  each  field  strength  were  13.3  ±  1.6  and  12.9  ±  2.3  s"'  mG’’  for  SPIO 
and  USPIO  respectively.  This  too  is  expected  if  cellular  bulk  magnetic  susceptibility  satisfies 
conditions  of  the  static  dephasing  regime  and  is  the  dominant  mechanism  of  relaxation. 

The  7?^*  relaxivities  of  both  SPIO  and  USPIO  suspended  in  gel  were  significantly  less  than 
those  of  the  iron-oxide  agents  compartmentalized  within  cells,  being  a  factor  of  3  and  8  less  than 
SD  regime  theory  predictions,  respectively.  The  reduction  in  7?^*  for  free,  uniformly  distributed 
SPM  particles  is  expected  since  individual  iron-oxide  cores  do  not  have  enough  iron  to  produce 
outer  sphere  inhomogeneity  patterns  which  are  strong  enough  to  prevent  the  rate-reducing  effects 
of  motional  averaging.  Iron-oxide  cores  have  a  mass  of  only  1.92  x  10’^  pg  of  iron.  This  is 
derived  using  the  known  density  of  iron-oxide  (3.76  gFe/cm^)  (28)  and  the  size  of  the  iron-oxide 
core  (4.6  nM)  (4,27).  In  the  appendix,  the  SD  regime  definition  is  reformulated  in  terms  of  LMD, 
in  mG,  and  iron  mass  per  compartment  («,.),  in  pico-grams.  The  result  is  Eq.  [A7],  shown  here: 
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For  single-core  iron-oxide  particles  in  uniform  suspension,  using  1.92  x  10"’  pg  as  the  mass  of 
iron  per  iron-oxide  core,  accurate  SD  regime  predictions  are  only  expected  for  LMD  »  10’^  mG. 
From  the  calibration  shown  in  Fig.  2,  this  corresponds  to  an  iron  concentration  greater  than  7  x 
10*  mg/ml,  which  is  much  larger  than  can  be  obtained  in  practice.  This  SD  regime  limit  describes 
the  case  of  SHU  555C,  the  USPIO  studied,  as  it  is  a  monocrystaline  iron-oxide  nanoparticle 
(MION),  having  only  one  iron-oxide  core  per  particle.  The  SPIO  studied,  SHU  555A,  is  a 
polycrystaline  iron-oxide  nanoparticle  (PION)  having  multiple  iron-oxide  cores  per  particle. 
Since  SPIO  has  a  greater  mass  of  iron  per  particle,  SPIO  in  uniform  suspension  has  a  larger  7?^* 
relaxivity  as  it  more  closely  approaches,  though  does  not  yet  reach,  the  SD  regime  limit  for 
practical  LMD  values,  as  expressed  by  Eq.  [4].  The  agreement  of  the  relaxivity  measurements 
for  the  cellular  suspensions  with  the  predictions  of  SD  regime  theory  is  striking,  though  upon 
reflection  not  unexpected,  since  table  1  indicates  that  the  smallest  mass  of  iron  within  any  cell 
used  was  0.7  pg.  In  this  case,  SD  regime  predictions  are  valid  for  LMD  »  0.1  mG.  All  iron- 
oxide  loaded  cell  samples  in  our  experiments  had  measured  LMD  values  well  above  this  range. 
We  conclude  from  these  experiments  the  remarkable  fact  that  cellularly  compartmentalized  iron- 
oxide,  at  any  reasonable  loading  level,  leads  to  Rf  relaxivity  that  is  both  maximized  and  can  be 
easily  and  accurately  predicted  from  SD  theory,  whereas  free  iron-oxide  nanoparticles  are 
neither  large  enough  nor  dense  enough  to  approach  the  SD  relaxivity  limit. 

To  further  investigate  the  reliability  of  static  dephasing  regime  predictions  for  cellularly 
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compartmentalized  iron-oxide,  the  CPMG  90-180  pulse  interval  (x^p^g)  dependence  of  was 
explored  at  4  T  for  the  constant  cellular  density  groups.  Figure  6  shows  the  7?^  versus  LMD 
curves  for  SPIO  and  USPIO  acquired  using  a  x<^„g  of  5  ms  and  10  ms.  As  detailed  in  the  Theory 
section,  no  explicit  analytic  theory  for  large  magnetic  perturbers  exists  to  predict  dependence, 
so  the  solid  lines  are  merely  second  order  polynomial  fits  to  the  data  to  guide  the  eye.  There  are 
two  clear  observations:  R2  increases  with  LMD  for  all  curves  and  R2  is  greater  for  the  x^p^g  =10 
ms  than  the  Xcp„,g  =  5  ms  curve.  There  are  two  possible  explanations  for  the  observed  increase  in 
R2  with  LMD\  intracellular  diffusion  through  outer  sphere  dipolar  field  inhomogeneities  of  iron- 
oxide  cores,  as  best  described  by  quantum  dipolar  relaxation  theory;  or  extracellular  diffusion 
through  dipolar  field  inhomogeneities  produced  by  cellular  BMS  differences.  The  second 
observation  of  increased  R2  with  x^p^g  indicates  that  BMS  effects  contribute,  at  least  in  part,  to  the 
R2  relaxation,  and  that  some  deviation  from  predictions  of  the  SD  regime  are  observable. 
However,  this  non-zero  BMS  contribution  to  R2  is  slight,  being  at  least  70  times  less  than  R2,  so 
deviation  from  the  purely  static  dephasing  regime  for  cellularly  compartmentalized  iron-oxide  is 
minimal. 

R2  sensitivity  to  compartmentalized  iron-oxide  concentration  was  fiirther  explored  by  varying 
concentration  through  changes  in  cell  density  versus  iron  content  per  cell.  Figure  7  plots  R2 
versus  LMD  for  the  variable  cell  loading  samples  (dose  and  incubation  time  groups)  and  the 
variable  cell  density  samples  (cell  density  group)  for  both  SPIO  and  USPIO  at  4  T  with  x^p^g  = 
Sms.  Again,  LMD  is  directly  proportional  to  iron  concentration  and  second  order  polynomial  fits 
are  included  merely  to  guide  the  eye.  is  clearly  more  sensitive  to  increased  iron  concentration 
for  the  variable  cell  density  group  than  for  the  variable  cell  loading  group.  Two  properties 
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explain  this  result:  sensitivity  is  a  measure  of  deviation  from  the  predictions  of  SD  regime 
theory;  and  SD  regime  predictions  become  more  accurate  as  the  mass  of  iron  per  cell  increases 
(see  Eq.  [4]).  For  the  variable  cell  loading  group,  increasing  iron  concentration  by  increasing  iron 
mass  per  cell  causes  the  samples  with  more  iron  per  cell  to  better  approximate  the  SD  regime  and 
hence,  to  have  reduced  relative  to  the  variable  cell  density  group.  This  argument  also  provides 
an  explanation  for  the  sub-linear  response  of  R2  to  LMD  for  the  variable  cell  loading  group. 

Discussion 

Cellular  Iron  Uptake 

Before  studying  the  relative  sensitivity  and  specificity  of  various  NMR  relaxation  parameters 
{R,,  R2  and  R2')  to  the  presence  of  iron-oxide  in  cells,  a  means  of  assessing  the  mass  of  iron 
ingested  through  phagocytosis  was  necessary.  The  RMM  method  of  MR  susceptometry  (24,25), 
modified  by  us  to  directly  image  the  field  inhomogeneity  pattern,  met  this  requirement.  In 
addition,  MR  susceptometry  served  to  measure  sample  magnetization  and  thus  allowed  a  direct 
assessment  of  SD  regime  theoretical  predictions.  Quantifying  iron-oxide  uptake  during 
phagoc5^osis  confirmed  that  the  magnetic  label  does  not  decay  over  the  48  hour  iron-oxide 
incubation  time  studied  (29).  Preservation  of  the  magnetic  label  over  a  number  of  days  was 
confirmed  by  the  agreement  between  AES  and  susceptometry  estimates  of  iron  mass  per  cell. 
Breakdown  of  the  iron-oxide  core  following  phagocytosis  would  have  destroyed  the  SPM  core 
magnetic  moment  causing  sample  susceptibility  to  drop  dramatically  (29).  Preservation  of  the 
SPM  magnetic  moment  after  phagocytosis  is  important  for  achieving  optimal  sensitivity  and 
quantifiability  of  macrophage  cell  number/density  in  vivo. 


19 


The  uptake  of  iron-oxide  was  found  to  be  linear  with  dose  and  incubation  time  to  very  large 
loading  levels,  permitting  the  definition  of  an  uptake  efficiency.  Linearity  with  incubation  time 
was  previously  demonstrated  by  Pratten  and  Lloyd  (30)  using  radioactively  labeled  polystyrene 
beads  with  peritoneal  macrophages.  Uptake  linearity  is  confirmed  in  this  work  for  dose  variation 
and  for  the  more  applicable  SPM  particles.  Uptake  of  iron-oxide  was  nearly  four  times  more 
efficient  with  the  larger  SPIO  particles  than  with  USPIO.  This  is  consistent  with  the  result  of 
Pratten  and  Lloyd  (30)  that  described  an  increased  phagocytic  contribution  to  particle  volume 
uptake  with  increasing  size  of  polystyrene  beads.  Since  the  iron  density  for  the  two  SPM 
particles  studied  differ,  a  quantitative  comparison  of  the  relative  uptake  rates  for  these  particle 
sizes  with  the  work  of  Pratten  and  Lloyd  is  not  possible.  However,  more  important  for  targeted 
imaging  applications,  the  mass  of  iron  ingested  within  cells  is  predictable,  with  a  unique  uptake 
efficiency  being  described  for  every  nanoparticle  and  phagocytic  cell  type. 

Being  able  to  quantify  relative  uptake  efficiency  for  different  SPM  particles  has  implications 
both  for  in  vitro  and  in  vivo  loading  methods.  In  vitro  loading  methods,  used  for  cell  trafficking 
(8),  require  fast,  efficient  loading  of  iron-oxide  into  cells  prior  to  re-injection.  Quantifying  uptake 
of  different  iron-oxide  particle  preparations  may  help  clarify  the  most  efficient  particle  for  this 
purpose.  Knowledge  of  the  iron-oxide  mass  per  cell  is  important  for  estimates  of  cellular 
detection  threshold.  In  vivo  loading  methods,  where  SPM  particles  are  injected  directly  into  the 
bloodstream,  may  also  benefit  from  quantification  of  iron-oxide  uptake.  There  are  many  factors 
affecting  the  efficiency  of  SPM  uptake  by  targeted  cells  (7).  These  include  RES  clearance  of 
SPM  particles  from  the  bloodstream,  tissue  perfusion,  target  cell  uptake  efficiency,  SPM  particle 
breakdown  and  SPM  particle  clearance  from  tissue  via  the  lymphatic  system.  Particle  clearance 
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and  tissue  perfusion  rates  can  be  measured  by  other  methods.  Therefore,  measurement  of  target 
cell  uptake  efficiency  may  provide  the  critical  missing  information  for  estimating  in  vivo  rates  of 
particle  uptake  by  targeted  cell  populations. 

Sensitivity  and  Quantification 

The  NMR  parameter  (R„  R2  or  R2')  which  was  most  sensitive  to  the  presence  of  iron-oxide 
within  cells  was  explored.  As  indicated  in  table  2, 7?^*  relaxivity  was  70  times  greater  than  7?^  and 
3100  times  greater  than  Rj,  for  iron-oxide  compartmentalized  within  cells.  The  main  effect  of 
compartmentalization  was  to  greatly  increase  the  relative  7?2*  sensitivity,  compared  to  the 
uniformly  distributed  iron-oxide  case,  where  7?^*  and  7?^  are  of  the  same  order  (4,14).  However, 
the  relative  sensitivity  of  image  signal  intensity  to  the  presence  of  iron-oxide  is  not  determined 
solely  by  relaxivity.  The  relative  contrast  to  noise  (CNR)  of  7),  T2  or  Tf  weighted  acquisitions, 
with  and  without  iron-oxide  present,  is  derived  from  iron-oxide-induced  relaxation  enhancement, 
as  well  as  the  baseline  tissue  relaxation  rates  measured  prior  to  iron-oxide  injection  (15).  While 
baseline  tissue  relaxation  rates  are  different  for  every  tissue  and  field  strength,  Rf  is  typically  no 
more  than  1  order  of  magnitude  greater  than  7?^  and  2  orders  of  magnitude  greater  than  7?y.  Since 
the  iron-oxide-induced  relaxation  enhancement  for  Rf  was  measured  to  be  approximately  2  and 
3  orders  of  magnitude  greater  than  for  7?^  and  R„  respectively,  7?^*  weighted  acquisitions  are 
potentially  the  most  sensitive  to  the  presence  of  cellularly  compartmentalized  iron-oxide 
particles.  Calculations  of  CNR  for  Rf  weighted  acquisitions  are  beyond  the  scope  of  this  work, 
but  are  possible  using  predictions  of  SD  regime  theory,  confirmed  for  SPM  loaded  cells  here, 
and  from  knowledge  of  baseline  tissue  72^*  sources. 
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A  new  approach  to  quantifying  the  concentration  of  iron-oxide  compartmentalized  within 
cells  is  another  outcome  of  this  work.  We  achieved  this  by  examining  which  exchange  regime 
equations  were  most  valid  for  describing  outer  sphere  dipolar  field  relaxation  in  the  case  of 
cellular  iron  compartmentalization.  The  quantum  solution,  which  relies  on  the  extreme  motional 
narrowing  approximation,  has  previously  been  shown  to  adequately  describe  the  Rj  and 
relaxivity  of  uniformly  distributed  SPM  particles  (12,13,31).  Since  this  theory  predicts  R2  and  R2* 
equivalence,  the  observation  that  R2*  is  as  much  as  70  times  larger  than  R2  relaxivity  for 
cellularly  compartmentalized  SPM  clearly  invalidates  predictions  of  the  quantum  solution  in  this 
case.  This  demonstration  of  homogeneous  broadening  on  the  NMR  timescale,  and  the 
dependence  of  R2,  are  consistent  with  cellular  scale  BMS  differences  being  the  mechanism  of  R2 
sensitivity  to  cellularly  compartmentalized  iron-oxide.  The  magnetic  moment  of  the  entire  cell 
taken  as  a  single  magnetic  compartment  produces  a  large  enough  outer  sphere  dipolar  field  to 
violate  the  extreme  motional  narrowing  assumption  implicit  in  the  quantum  theory. 

Though  much  smaller  than  i?;*,  and  without  a  clear  theoretical  framework  to  predict 
relaxation  rate  for  cellularly  compartmentalized  iron-oxide,  R2  does  show  sensitivity  to  iron 
concentration  and  so  is  a  candidate  for  empirical  quantification  of  iron  concentration.  As 
indicated  in  Fig.  7,  R2  demonstrates  iron-oxide  sensitivity,  which  is  different  depending  on 
whether  iron  mass  per  cell  or  cell  concentration  is  varied  to  alter  sample  iron-oxide 
concentration.  Concentrating  more  SPM  into  fewer  cells  produces  a  better  approximation  to  the 
SD  regime  and  hence,  a  reduction  in  R2.  Therefore,  R2  is  sensitive  to  both  concentration  and 
distribution  of  SPM,  which  greatly  complicates  its  use  for  quantifying  iron-oxide  concentration. 
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The  SD  regime  predictions  for  R2'  agreed,  within  uncertainty,  with  measurements  for  the 
cellularly  compartmentalized  SPM  in  two  of  four  cases  (SPIO  at  4  T  and  USPIO  at  1.89  T), 
while  measured  was  significantly  higher  than  SD  regime  theory  for  the  other  cases  (USPIO  at 
4  T  and  SPIO  at  1 .89  T).  Overall,  we  observe  a  21%  elevation  in  relative  to  SD  regime  theory 
for  a  volume  fraction  of  0.189%.  Brooks  et  al.  (31)  have  noted  a  33%  increase  in  R2*  with  Monte 
Carlo  simulations  compared  to  SD  regime  theory  predictions.  The  authors  postulate  fliat  the  large 
perturber  volume  fraction  in  their  simulations  (0.5%)  violated  the  dilute  perturber  approximation 
implicit  in  SD  regime  theory.  In  this  case,  field  profile  overlap  between  different  perturbers 
eliminates  signal  from  regions  furthest  from  the  perturber,  where  field  gradients  are  the  lowest, 
leading  to  an  increased  relaxation  rate.  Monte  Carlo  simulations  with  very  small  volume  fraction 
(2  X  10  ®),  by  Muller  et  al.  (17),  produce  R2'  values  that  are  elevated  by  10%,  relative  to 
theoretical  values,  but  which  agree  within  uncertainty  (as  noted  by  Yablonskiy  and  Haacke  (22)). 
For  most  in  vivo  applications,  the  volume  fraction  of  cells  containing  iron-oxide  will  likely  be 
small  enough  to  satisfy  the  dilute  perturber  requirement;  hence,  we  expect  that  SD  regime  theory 
will  be  more  closely  approximated  in  vivo  than  in  our  more  artificial  cell  culture  conditions. 

As  detailed  in  the  results  section,  all  iron-oxide  loaded  cell  samples  tested  met  the  SD  regime 
requirement  for  mass  of  iron  per  compartment,  while  the  uniformly  distributed  SPM  in  gel 
samples  did  not.  As  a  result,  uniformly  distributed  SPM  samples  had  R2’  relaxivities  3  and  8 
times  less  than  SD  regime  theory  for  SPIO  and  USPIO,  respectively.  Although  samples  which 
clearly  satisfied  or  clearly  failed  the  static  dephasing  regime  condition  were  measured  in  this 
work,  the  transition  region  defining  onset  of  the  SD  regime,  given  in  Eq.  [2],  was  not  carefully 
probed.  However,  Monte  Carlo  simulations  by  Muller  et  al.  (17)  do  support  the  SD  regime 
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equation  defined  by  Yablonskiy  and  Haacke  (22).  Equation  [2]  defines  the  SD  regime  as  > 
0.19  |xm,  while  the  conventional  definition,  given  as  Rf  »  6Z)/5©s,  defines  the  SD  regime  as 
>  0.02  |im  (calculations  using  the  parameters:  f^  =  2x  10'*,  5®s  =  3.4  x  10’  s‘'  and  D  =  2.5x  10'® 
cm’  s'*).  Figure  1  in  the  work  of  Muller  et  al.  supports  the  prediction  of  Eq.  [2]  since  has 
clearly  reached  the  peak  value  for  >  0.19  jim,  while  is  only  60%  of  peak  for  R„  >  0.02  p,m 
in  that  work. 

For  cellularly  compartmentalized  SPM  particles,  R2'  is  potentially  the  most  sensitive  and  the 
most  quantifiable  NMR  parameter.  However,  it  is  necessary  that  the  mass  of  iron  per  cell  be 
sufficient  to  guarantee  static  dephasing  regime  predictions  in  order  to  relate  R2'  with  iron-oxide 
concentration.  Assuming,  for  the  sake  of  argument,  that  the  presence  of  SPM  particles  must 
enhance  R2'  by  at  least  100  s'*  to  be  well-distinguished  from  background  sources,  LMD  and  [Fe] 
must  exceed  9.3  mG  and  6.9  fxg/ml,  respectively.  This  is  derived  from  Eq.  [1],  and  the 
calibration  curve  in  Fig.  2,  respectively.  For  this  minimally  detectable  R2,  the  mass  of  iron  per 
particle,  as  calculated  using  Eq.  [4],  must  exceed  0.07  pg  per  cell  to  guarantee  SD  regime 
predictions.  Based  on  loading  levels  that  we  and  others  have  achieved  (4),  we  expect  that  this 
condition  will  be  met  by  most  cell  trafficking  studies  (which  use  in  vitro  loaded  cells)  and  may 
be  met  for  studies  of  inflammatory  models  using  injection  of  free  SPM  particles  with  subsequent 
in  vivo  cell  uptake. 

Conclusions 

The  relaxation  properties  of  SPM  particles  that  are  compartmentalized  within  cells  were 
studied  and  found  to  satisfy  predictions  of  static  dephasing  regime  theory  under  conditions 
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typical  for  many  in  vivo  applications.  The  required  conditions  are  that  the  mass  of  iron  per  cell 
exceed  a  known  threshold  to  ensure  the  validity  of  SD  regime  predictions.  To  test  SD  regime 
theory,  it  was  necessary  to  develop  methods  measuring  cellular  iron-oxide  uptake  and  sample 
susceptibility.  An  imaging  version  of  the  RMM  susceptometry  technique  was  developed  to  meet 
these  requirements.  In  addition,  the  linearity  of  phagocytic  uptake  of  iron-oxide  with  dose  and 
incubation  time  was  demonstrated  for  THP-1  cells,  permitting  the  definition  of  an  uptake 
efficiency,  unique  for  every  cell  and  particle  type. 

Under  conditions  satisfying  the  SD  regime,  the  greatest  sensitivity  to  cellularly 
compartmentalized  iron-oxide  was  observed  with  rather  than  R2  or  Rj.  This  differs  from  the 
case  of  uniformly  distributed  SPM  particles,  where  relative  sensitivity  between  parameters  is 
more  similar.  In  addition,  R2'  was  best  for  quantifying  iron-oxide  concentration.  This  is  because 
analytical  theories  for  and  R,  are  lacking  when  SPM  particles  are  compartmentalized  within 
cells,  whereas  SD  regime  theory  correlating  SPM-particle-induced  R2  enhancement  with  iron- 
oxide  concentration  was  proven  to  be  accurate.  R2'  showed  sensitivity  to  iron-oxide  concentration 
alone,  rather  than  both  concentration  and  mass  of  iron  per  cell,  as  was  observed  with  R2.  The 
non-unique  dependence  of  R2  on  mass  of  iron  per  cell  greatly  complicates  its  use  for  quantifying 
cellularly  compartmentalized  iron-oxide  concentration.  For  in  vivo  conditions,  R2' may  yet  prove 
to  be  both  the  most  sensitive  and  the  most  quantifiable  parameter  for  estimating  iron-oxide 
concentration. 

Appendix 

The  static  dephasing  regime  definition,  formulated  by  Yablonskiy  and  Haacke  (Eq.  [2]),  can 
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be  rewritten  as: 
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where  definitions  listed  in  the  introduction  for  ScOj  and  LMD  (LMD^^J,  in  SI 

units,  have  been  used.  Solving  for  f^,,  one  obtains  the  expression: 
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The  inequality  involving  f„  has  reversed  from  the  initial  expression  reported  in  Eq.  [Al].  This 
form  is  more  intuitive  since,  for  a  given  local  magnetic  dose  (iron  concentration), 
compartmentalization  of  magnetic  particles  to  less  than  some  maximum  fractional  volume,  f^,  is 
necessary  for  accurate  static  dephasing  regime  predictions. 


An  expression  that  describes  the  static  dephasing  regime  in  terms  of  the  required  mass  of  iron 
per  compartment  is  desired.  An  expression  calibrating  magnetic  dose,  LMD,  with  iron 
concentration  is  necessary  to  achieve  this: 

LMD  =  [A3] 


where  n^  is  the  compartment  density,  m^  is  the  mass  of  iron  within  each  compartment  and  is  a 
calibration  coefficient  relating  iron  concentration  to  magnetic  dose.  Using  the  fractional  volume 
definition,  fySn^V^  where  is  the  compartment  volume,  and  Eq.  [A3],  Eq.  [A2]  becomes: 
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Canceling  common  powers  of  n^  and  and  squaring  both  sides,  the  equation  simplifies  to: 
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where  LMD  has  been  reinserted  using  Eq.  [A3].  Finally,  the  static  dephasing  regime  equation 
can  be  evaluated  for  either  m,,  or  LMD  as  follows: 
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The  SD  regime  predictions  will  only  be  reliable  if  the  mass  of  iron  per  compartment  is  great 
enough  to  satisfy  Eq.  [A6]  for  the  minimum  LMD  in  the  range  being  measured.  Conversely,  for 
a  given  m^,  measurements  are  reliable  only  when  LMD  exceeds  that  in  Eq.  [A6]. 

The  expressions  in  Eq.  [A6]  can  be  further  simplified  if  the  calibration  coefficient  between 
magnetic  dose  and  iron  concentration,  Kpe,  is  known.  In  Fig.  2,  this  calibration  was  shown  to 
agree  within  uncertainty  for  the  two  iron-oxide  particles  studied,  with  a  mean  value  of  Kp^  = 
1.340  X  10'^  mG/(pg/mm^).  Using  this  value  for  Kp^,  the  self-diffusion  coefficient  for  free  water 
(D  =  2.5  X  10’^  mm^  s"')  and  the  proton  gyromagnetic  ratio  (y=  26.752  rad  mG‘' s"'),  Eq.  [A6]  can 
be  simplified  as  follows: 
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where  the  mass  of  iron  within  a  compartment,  m„  is  in  pico-grams  and  the  local  magnetic  dose, 
LMD,  is  in  milli-gauss. 
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Table  1:  Summary  of  Cellular  Iron  Uptake  Data 


SPIO  in  Cells 

USPIO  in  Cells 

LMD- 

[Fojsampifl  ** 

mass/cell " 

LMD- 

[F6]sample  ** 

mass/cell  “ 

(mG) 

(ug/ml) 

(pg  of  Fe) 

(mG) 

_ (ug/ml) _ 

(pg  of  Fe) 

(mg/ml) 

Dose  Group  (24  hours,  4  million  cells/ml) 

0.05 

31.8 

24.1 

6.0 

7.4 

5.4 

1.4 

0.1 

46.4 

35.2 

8.8 

15.5 

11.4 

2.9 

0.2 

99.5 

75.4 

18.9 

23.9 

17.6 

4.4 

0.3 

129.2 

98.0 

24.5 

33.4 

24.6 

6.1 

(hours) 

Incubation  Time  Group  (0.1  mg/ml,  4  million  cells/ml) 

6 

19.1 

14.5 

3.6 

4.1 

3.0 

0.7 

12 

35.8 

27.2 

6.8 

6.7 

4.9 

1.2 

24 

46.4 

35.2 

8.8 

15.5 

11.4 

2.9 

36 

65.6 

49.7 

12.4 

18.2 

13.4 

3.3 

48 

76.2 

57.7 

14.4 

22.7 

16.7 

4.2 

(million  cells/ml) 

Ceil  Density  Group  (0.1  mg/ml,  24  hours) 

1 

13.3 

10.1  {9.5}- 

10.1 

6.8 

5.0  {4.0}- 

5.0 

2 

24.5 

18.6 

9.3 

10.0 

7.4 

3.7 

4 

46.4 

35.2  {46.4} 

8.8 

18.3 

13.5  {16.7} 

3.4 

8 

88.4 

67.0 

8.4 

28.3 

20.8 

2.6 

16 

167.0 

126.7  {170.9} 

7.9 

69.5 

51.1  {52.7} 

3.2 

“Local  Magnetic  Dose  of  cellular  samples  obtained  using  susceptometry  at  4  T.  Reported 
LMD  is  given  relative  to  the  reference  cellular  sample  that  was  not  exposed  to  iron-oxide. 
‘’Sample  iron  concentration  determined  from  LMD  and  calibration  in  Fig.  2. 

'Sample  iron  concentration  determined  from  AES  in  brackets. 

‘‘Mass  of  iron  per  cell  determined  from  sample  iron  concentration  and  mean  cell  count. 
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Table  2:  Relative  Relaxivities  at  1.89  T  of  Iron-Oxide  Loaded  Cells 


Relaxivities 

SPIO 

USPIO 

at  1.89  T 

(s  '  mG  ■') 

(s  ■'  mG  ■') 

Ri/LMD 

0.0036  ±0.0005  0.0051  ±0.0015 

Rz/LMD 

0.218±0.018 

0.165  ±0.023 

RzVLMD 

15.6±1.0 

11.3±2.1 

Table  3:  Effect  of  Compartmentalization  on  R2*  Relaxivity 


Compartment 

Ra'atl.SQT 

R2*at4T 

Type 

(s"'  mG‘') 

(s-'  mG  -’) 

SD  Theory  ® 

10.78 

10.78 

SPIO  in  Cells 

15.6±1.0 

11.0  ±1.3 

USPIO  in  Cells 

11.3±2.1 

14.5±1.0 

SPIO  in  Gel 

— 

3.8  ±0.3 

USPIO  in  Gel 

— 

1.4±0.2 

'Outer  Sphere  Static  Dephasing  Regime  prediction  from  Eq.  [1]. 


Figure  Captions: 

Figure  1:  Demonstration  of  susceptometry  method,  (a)  magnitude  image  and  (b)  field  map 
obtained  from  axial  slice  through  NMR  tube  concentrically  located  within  water  filled  graduated 
cylinder.  4  radial  field  map  profiles  were  selected  (heavy  bar),  plotted  individually  (c),  averaged 
together  (triangles  in  d),  and  fit  to  Eq.  [3]  (solid  line  in  d)  to  obtain  each  sample  susceptibility 
(Xi)-  The  estimate  for  Xi  was  -7.61  ±  0.02  ppm  for  the  0.05  mg/ml  SPIO  suspension  shown. 

Figure  2:  Calibration  curves  relating  SPM  magnetization  to  iron  concentration  for  SPIO 
(triangles)  and  USPIO  (circles)  suspensions.  The  LMD  is  defined  relative  to  pure  gelatin.  Linear 
regression  of  the  SPIO  (solid  line)  and  USPIO  (dashed  line)  data  produced  calibration  slopes  of 
1319  ±  69  and  1360  ±  55  mG  per  mg/ml,  respectively. 

Figure  3:  Cellular  iron  uptake  curves  indicating  the  mass  of  iron  per  cell  which  is  taken  up  per 
unit  dose  and  incubation  time.  Data  for  SPIO  (triangles)  and  USPIO  (circles)  are  shown  for  both 
the  dose  groups  (solid  symbols)  and  incubation  time  groups  (open  symbols).  The  slope  from 
linear  regression  of  the  SPIO  (solid  line)  and  USPIO  (dashed  line)  data  defined  the  uptake 
efficiency  as  3.06  ±  0.40  and  0.80  ±  0.08  pg  of  iron  per  mg/ml»hrs  respectively. 

Figure  4:  Relaxivity  curves  indicating  the  R,  (squares),  R2  (circles)  and  R2*  (triangles)  relaxation 
rates  as  a  function  of  LMD  for  cellularly  compartmentalized  SPM  samples.  Data  for  SPIO  (solid 
symbols)  and  USPIO  (open  symbols)  are  shown  for  the  combined  dose  and  incubation  time 
groups.  The  slope  from  linear  regression  of  the  SPIO  (solid  line)  and  USPIO  (dashed  line)  data 
defined  the  relaxivities  as  listed  in  table  2.  Note  that  R2*  is  nearly  2  orders  of  magnitude  greater 


34 


than  i?2  for  cellularly  compartmentalized  SPM. 


Figure  5:  R2*  relaxivity  curves  at  4  T  demonstrating  the  effect  of  SPM  compartmentalization.  The 
R2*  relaxation  rate  for  SPIO  (triangles)  and  USPIO  (circles)  are  shown  as  a  function  of  LMD  for 
cellularly  compartmentalized  SPM  samples  (solid  symbols)  and  free  SPM  suspensions  in  gel 
(open  symbols).  The  slope  from  linear  regression  for  the  SPIO  (thin  solid  lines)  and  USPIO  (thin 
dashed  lines)  data  determined  the  relaxivities  as  listed  in  table  3.  Note  that  R2*  increases 
dramatically  when  SPM  is  compartmentalized  within  cells  and  agrees  well  with  predictions  of 
SD  regime  theory  (thick  solid  line). 

Figure  6;  Demonstration  of  diffusion  effects  on  R2  relaxation  rate  sensitivity.  The  R2  relaxation 
rate  for  cellularly  compartmentalized  SPIO  (triangles)  and  USPIO  (circles)  are  shown  as  a 
function  oi  LMD  for  CPMG  acquisitions  having  T^pmg  equal  to  5  ms  (solid  symbols)  or  10  ms 
(open  symbols).  The  curves  for  SPIO  (solid  line)  and  USPIO  (dashed  line)  have  no  theoretical 
basis  and  serve  only  to  guide  the  eye.  Note  that  R2  increases  slightly  with  x^pnig  indicating  some 
deviation  from  the  static  dephasing  condition. 

Figure  7:  R2  relaxation  rate  curves  demonstrating  the  effect  of  SPM  compartmentalization.  The 
R2  relaxation  rate  for  SPIO  (triangles)  and  USPIO  (circles)  are  shown  as  a  function  of  LMD  for 
samples  having  constant  iron  mass  per  cell  but  variable  cell  density  (solid  s)mibols)  as  well  as 
those  having  constant  cell  density  but  variable  iron  mass  per  cell  (open  symbols).  The  curves  for 
SPIO  (solid  line)  and  USPIO  (dashed  line)  have  no  theoretical  basis  and  serve  only  to  guide  the 
eye.  Note  that  R2  sensitivity  to  SPM  is  much  larger  for  the  variable  cell  density  then  the  variable 
cell  loading  group  despite  having  the  same  sample  iron  concentration  for  each  LMD  value. 
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